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The spectral functions of t J and t Jxy models in the limit 
of J/t — * and at finite temperatures T <C t are calculated us- 
ing the spin-charge factorized wave function. We find that the 
Luttingcr-liquid like scaling behavior for a finite system with 
L sites is restricted below temperatures of the order T < J/L. 
We also observe weight redistribution in the photoemission 
spectral function in the energy range t, which is much larger 
than the temperature. 

Single-particle spectral functions are very useful to un- 
derstand the electronic structure of solids. They are mea- 
sured in photoemission [B(k, u>)] and inverse photoemis- 
sion h4(fc,u>)] experiments. For actual calculations, the 
Lehmann representation is very useful: 



B(k,u) = ±^\(f\a 

i>f t o- 



k,<7 1 



\i) S(lu - Ei+ E f )e 



(1) 



where i and / denote the initial and final states with N 
and N — 1 electrons, respectively, and a^.a annihilates 
an electron with momentum k and spin a. Furthermore, 
Z = J2i e ~^ Ei i s the partition function with (3 = 1/T 
being the inverse temperature. A similar expression holds 
for A{k,uo), which we will not treat in this paper. 

In contrast to quasiparticles in usual three-dimensional 
Fermi liquids, the collective^, excitations of onc- 
dimensional interacting electrons^ give rise to anomalous, 
scaling behavior of the one-particle Green's functional 
with non-universal exponents. For example, the momen- 
tum distribution function — ^ J du>B a (k,ui) takes the 
form nk ~ Tip + sgn(fc — kp)\k — kp\ a near the Fermi mo- 
mentum kp and zero temperature, where the exponent 
a depends on the actual model and coupling constants. 
Similarly, the local spectral function (single-particle den- 
sity of states) B(u>) = ■j^^ZkB{k,ijj) also scales with a 
power law:B B{u>) oc \w — ep\ a . To describe this critical 
behavior of one-dimensional models at low energies, Ha' 
dane introduced the fruitful concept of Luttinger liquids. 
Following a different approach, conformal field theory 
tells us that the exponents, are related to the finite size 
corrections of the energy BeI 

Recent experiments on quasi one dimensional matesi 
als raised the question if this behavior can be observedB 
Furthermore, in these experiments an anomalous spectral 
weight transfer has been observed: changing the temper- 
ature by 100 K, one can observe weight redistribution on 
the scale of 1 eV, which is a hundred times larger than 
the temperature. In this paper we will try to explain this 
behavior in a simple way. 

We are considering the isotropic and anisotropic t J 
model, defined by the Hamiltonian 



ff« = -*I>l )( A+i )CT + H.c.) 

+E E ja ( sfSi+i - 3<w^+0 , 

i a—x,y,z 

in the limit of small exchange J a — ► 0, where (x% (j are the 
usual projected operators to exclude double occupancy. 
Actually, the Hubbard model in the large-?/ limit can be 
mapped onto a strong coupling model usually identified 
as the tJ model plus three-site terms using a canonical 
transformation,!^ where J = At 2 /U is small. The spec- 
tral function of the Hubhaxd model has been studied uSi. 
ing exact diagonalizationlij and Quantum Monte Carloliil 
techniques, which both have well known limitations. 

An alternative, powerful but model limited approach is 
based on the special property of the wave functions of the. 
Hubbard model in the limit of large Coulomb repulsionllHl 
(also for J/t — > in the t J model), that the wave function 
factorizes: 



\i) = W) ® 1x0. 



(2) 



This has allowed the calculation and confirmation of the 
power law behavior of the static correlation function 
rik and gaveEj a = 1/8. \ip f ) describes the charge de- 
grees of freedom and is a wave function of free spinlcss 
fermions with momenta kp quantized as Lk'^ = 27rXj+Q', 
where Xj are distinct integers, < < L — 1 and 
j = 1,2,..., N. Twisted boundary conditions are im- 
posed by the momentum Q' of the spin wave function 
\x'), which describes the spins on a squeezed lattice of N 
sites and are eigenfunctions of a spin Hamiltonian with 
an effective spin exchange J which depends on the actual 
charge wave function \tp'), and e.g. for the ground state 
J a = J a n[l - sin 2 (7m)/(7rn) 2 ], where n = N/L. We will 
take periodic boundary conditions to avoid edge effectsE-a 
and an even number of electrons not a multiple of four 
(i.e. N = 2,6, 10, . . .) for convenience. 

To calculate the thermal average, we need to know all 
the energies and wave functions of the spin model. Since 
for the Heisenberg model this is very difficult to obtain, 
we turn to the XY model (i.e. J z = 0). In this spe- 
cial case, the spin model can be mapped onto noninter- 
acting spinless fermions using the Wigner-Jordan trans- 
formation. Assuming that the occupied sites represent 
the T spins, the states are characterized by Nf integer 
numbers > J'^ > N — 1, and the momenta q'j of the 
free spinless fermions representing the spins are quan- 
tized as Lq'j = 2'nJ'y Finally, the momentum of the spin 
wave function determining the boundary condition of the 
charge part is Q' — Yljli 1j = 2-nJ'jN, with J' integer. 
The energy of the state is simply = _E; iC + £j )S , where 
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E itC = ~2tJ2f =1 cosk' 3 and E itS = JxYj2fli C0S 1j, 



while the momentum reads Pi = Y] - —1 fc^. One should 
note that despite the fact that both the charge and 
spin wave functions in Eq. (|^) are those of free spinless 
fermions, the resulting wave function describes a nontriv- 
ial and strongly correlated system. As far as the exponent 
a (at T = 0) is concerned, it changes from a_j= 1/8 in 
the isotropic case to a = 1/4 in the XY casc.Ej 

Similarly, the final, N — 1 electron wave function fac- 
torizes as well: |/) = \ifj) ® \x)- The quantum numbers 
for the spinless fermions representing the charges are Ij , 
and the corresponding momenta Lkj = 2nTj + Q. Here 
Q = 2-kJ/{N — 1) is the momentum of the N — 1 spin 
wave function, (0 > J > N - 2). 

Since the charge and the spin part are coupled through 
the momentum Q' of the spixLwave function, the partition 
function does not factorizeEj (i.e. the free energy is not 
a sum of charge and spin contribution) and it will read 
Z = Eq>Zc(Q')Z s (Q'), where 



Z C (Q') 



E« 

{!'} 



-0E i: 



z s (Q') 



E 



-PEi,< 



and the sum in Z s is over the states with given momen- 
tum Q' . In calculating the thermodynamic averages, one 
can work in principle in an ensemble fixing either the 
magnetization or the magnetic field. We have used both 
ensembles, and although the results in the thermody- 
namic limit should be independent of the ensemble we 
choose, there are strong finite size effects. 

Even though we know all the excitations for the 
tJxY model, we will make further restrictions which are 
needed to perform calculations on reasonably large sys- 
tem sizes: Namely, we will consider temperatures much 
smaller than the energy scale of the charges. In other 
words, for the charge part we neglect the excitations 
and take the ground state given by consecutive inte- 
gers {!'} = {-N/2, -N/2 + 1, . . . , N/2 - 1}. Then the 
remaining free parameter is T / J, and all the tempera- 
ture dependence is now in the spin part. Furthermore, 
since the energy of the charge part also depends on Q' as 
Ei. c {Q')-E^ c (Q' = 7r) = 2^u c (Q'-7r) 2 , where u c cx t is 
the charge velocity, we will assume that the momentum 
of the spin part in the initial N electron state is Q' = tt. 
This restriction is actually more for convenience, as the 
result does not depend on this assumption - we will com- 
ment on this later on. 

Using the factorized wave function, the spectral func- 
tion defined in Eq. ([!]) simplifies tcO 



(3) 



Q,(T 



Here Bq{k,uj) — _BQ Q/ =7r (A;, to) depends on the spinless 
fermion wave function only: 

Bq,Q>{k,uj) = L^2\{ip Q \b \ipQ, 
{1} 

x5(u- E tiC + Ef >c )5k,p,-Pf, 



where bo annihilates a spinless fermion at site 0. The 
matrix elements in BQ.Qi(k,u) read: 

T 1/ / |i \ ij \| 2 7--2N+2 • 2N-2 Q ~Q 

L\(ipQ\b \ipQ,)\ =L sm 



IT 



x I I sm 

j>i 



2 kj ' 



- kj 



j>i 



V - V 



2 
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We can actually recognize Anderson's orthogonality 
catastropheEj in these complicated matrix elements, 
which is a consequence of changing the boundary con- 
dition from Q to Q' in the charge wave function due to 
momentum transferred to the spins. 

On the otheij hand, the contribution of the spin degrees 
of freedorrMiJ D a (Q, (3) = D a (Q, Q' = tt, (3) is given by 



D a {Q,Q',p) 



1 



1 



N-1Z„ 



E 

,{J'}q> 



im(Q'-Q) -/3Ei 



where wo^m.o- denotes the amplitude to transfer a a spin 
from site to m: 



i,<7 = (x'\Pm, m -i ■ ■ ■ Pi,on \x')- 



(4) 



The operator Pj.j+i = 2SjSj + i + | permutes the spins 
on sites j and j + 1. 

a. Spin part: For the XY model, after introducing 
the spinless fermions (with operators /) in the Wigner- 
Jordan transformation, the permutation operator reads 



p j+i,] = n 3+i n j + fj+Jj + (i-nj+iXl-nj) + f]f j+1 , 

and the spin transfer amplitude can be easily calculated 
from Eq. (0) using Wick's theorem. We find that 
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where gi = (-l)'(x|///o|x) = £ Efti e^'^K In par- 



N 

ticular, go — N^/N and g-j = g*, furthermore the rela- 
tion U} -tN-l-m,<T = e tQ 'uo->m,a holds. 

For large temperatures t ^> T 3> J (equivalent to 
"hot spins" of Ref. |2^) a high temperature expansion is 
possible: if we relax the constraint that we take only 
states with momenta Q' , then it follows that Z s = 
2 N + 0((3Jxy) and for LUo-* m ,a we have to count the 
number of states where the first m + 1 spins have S z =f , 
which is 2 N - m - 1 . Working in a subspace with definite 
momentum Q' , each subspace will acquire roughly l/N 
of the values given above (the actual distribution depends 
on how many states are in a given Q' subspace), and in 
the thermodynamic limit we get 



D a (Q,Q',f]^0) 



1 



N-l 10-8cos(Q-Q') 



(5) 
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This result is not only valid for the XY model, but also 
for the isotropic Heisenberg model. 

We show the behavior of the spin part in Fig. El. Apart 
from the clear power-law singularity near Q = ir/2 at 
zero temperature, we observe that at fixed small tem- 
perature this behavior disappears as we increase the sys- 
tem size. This indicates that the singularity will vanish 
for any finite temperature in the thermodynamic limit. 
Also there is a difference between calculating D a (Q,f3) 
in the two ensembles mentioned above, however the finite 
size effects are decreasing with increasing N. Let us also 
note that the sum rule D a {Q, (5) = N a /N is satis- 
fied for any temperature. Furthermore, D a (Q,Q' , (3) « 
D a {Q — Q' + 7r, 0) in the thermodynamic limit. 




X AT - 1 V r i(Q-Q')(N l -,n) 

ONi-m - N _ 1 2^ e 

Q 

and comparing Eqs. (^) and (|t]), we get 

B Q .Q>{k) =J2 em W\4 II ^ nAQ - Q,) c Q W). (8) 
I l'=\ 

This can be further simplified using the identity 
e mi(Q-Q ) _ 2 + where we introduced a; = e 1 ^ - ^ ) — 

1, so that (ip'\c] U\~2i e m ' ,(Q ~ Q ' ) c |V') in Eq. (§ is equal 
to 

(4 c o) ( c / c i) ■■■ (4 c i-i) 

x{c\c ) l + x(c\ Cl ) ... xiclc^) 
x(c\c Q ) x(4c t ) ... xiclc^) 



x ( c i-i c o) x ( c i-i c i) ■■■ l + x ( c l-l c l-l) 

where {c\c v ) = J^ . e~ tk ^ l ~ l ). Using this equation, 
we are able to compute BQ(k 1 Lu) for systems with a 
few hundred sites. It turns out that BQ^Qi(k,oj) = 
i?Q_Q/ +7r (fc, ui) apart from some small finite size correc- 
tions, therefore our assumption to fix Q' = 7r is justified. 
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FIG. 1. Temperature dependence of D a (Q,/3) for the XY 
model in zero magnetic field (solid) and zero magnetization 
(empty symbols) for T/J = 0, 0.5 and T > J. The solid line 
for T = shows the N = 250 result, and for T > J Eq. (|) 
is plotted. 

b. Momentum distribution function: From Eq. (||) we 
get: 

n k = Y,B Q (k)D a (Q,f3), (6) 
Q 

To calculate efficiently, we have to find a convenient 
way to evaluate Bq(/c). For that reason, let us follow 
Ref. [l3]: In the alternative representation of the momen- 
tum distribution 

i 1=0 

we replace \i) by the factorized wave function, Eq. (^|): 

JV-2 

(*l o i,a ffl 0, ff l i > = ^2(lp'\4 5 fft-mC \lp'Wo-,m,<T, (7) 
m=0 

where N[ = J^;'=o n /' coun t s the number of spinless 
fcrmions between sites and I, and loq *m,& is calcu- 
lated for the particular \x')- Now, replacing 5jv t _ m by 
its Fourier representation: 




7t/4 71/2 371/4 7C 

k 

FIG. 2. Momentum distribution of the tJxY model for 
T = (solid line: L = 500), T/J = 0.5, and ( » T » j 
(solid line: L = 300) for quarter filling (L — 2N) and zero 
magnetic field. There are strong finite size effects for T — J/2. 
In the insert we show the scaling of Anj? for zero magnetic 
field (solid) and zero magnetization (empty symbols) obtained 
from different system sizes and temperatures. 

We show our numerical results in Fig. |[ The T = re- 
sult shows power law behavior at the Fermi momentum. 
Increasing the temperature, the power law survives until 
some crossover temperature (~ J/L), where it becomes 
a continuous function of momentum for higher temper- 
atures. To study this behavior in detail, we concentrate 
on the jump at kp, defined as Anp = n k - — n k +, where 

kp = kp ± ttN/L are the momenta of the finite system 
closest to the Fermi point. This jump is finite for finite 
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size system and scales with L~ a in the Luttinger-liquid. 
If the singularity disappears and n}~ becomes a continu- 
ous function around then An F oc 1/L. In the inset of 
Fig. ||we show the "size independent" L a An F vs. LT/J. 
It is remarkable, that at low temperature the points fol- 
low a universal curve: 

An F = L~ a f(LT / J). 

A crossover temperature, scaling with J/L, can be clearly 
observed, and for larger temperatures L a Anp — > 0. This 
behavior can be understood if we recall that the tem- 
perature enters by dividing the energy oc 2iru a /L of the 
low-energy excitations. 

c. Local spectral function: The single-particle density 
of states is given by 



(9) 



where Bq{uo) — j- Bq{k, uS). Let us concentrate on 
the isotropic t J model (equivalent to large-?/ Hubbard 
model) in the limiting T = and t ^> T 3> J cases only. 
At low temperatures D a (Q,(3) is large near Q = tt/2, 
and the largest part in the convolution (^|) comes from 
Bg =7r / 2 (w). For the "hot spin" case, D a {Q,f3) is large 
near Q = tt, and Bq = ^{uj) gives most of the contribution 
to B(u>), shown in Fig. [|. In other words, increasing the 
temperature we transfer less and less momentum to the 
spins, and the role of the orthogonality catastrophe in 
Bq{u>) decreases. Since changing Q results in a consider- 
able redistribution of the weight in Bq(lu) (see the inset 
in Fig. |J), the weight transfer of B{uo) at the energy scale 
of t is due to the temperature dependence of D a (Q,[3) 
set on a much smaller temperature scale - naively we 
would expect smearing of B{w) near the Fermi energy 
within \uj — e F \ « T. We should also note that the di- 
vergence of the spectral function at the Fecmi energy is 
purely the artifact of the J/t — » limittll. For finite 
J, the local spectral function has a broad peak around 
u> s=s 77 J due to the spinon dispersion, and a second broad 
peak near the band edge (ui ~ 2t). The weight transfer 
then would be from the "spinon" to the "holon" peak. 
A similar weight redistribution, is observed in the two- 
dimensional t J model as well.EJ 
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FIG. 3. Local spectral function for T — (solid) and 
t S> T ^> J (dashed line) for the quarter-filled Hubbard 
model, L — 220. For this particular filling ef — 0. In the 
inset: Bq(u>) for different values of Q. 

To conclude, we have studied the temperature evolu- 
tion of the momentum distribution function and local 
spectral function. First, we give a method to calculate 
rife for large system sizes for the tJxy model at zero tem- 
perature. Next, we observed that the power-law behavior 
is restricted to temperatures inversely proportional to the 
system size. In the thermodynamic limit the system is 
critical at T = only. Finally, a weight redistribution 
in the single-particle density of states takes place over 
a broad energy range, which can be easily understood 
using the concept of "spin-charge" separation. 

We would like to thank H. Frahm, J. Jaklic, H. Shiba 
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